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Abstract 

Results are presented of a full leading-order evaluation of the shear viscosity, flavor diffusion con- 
stants, and electrical conductivity in high temperature QCD and QED. The presence of Coulomb 
logarithms associated with gauge interactions imply that the leading-order results for transport 
coefficients may themselves be expanded in an infinite series in powers of l/log(l/g r ); the utility of 
this expansion is also examined. A next-to-leading-log approximation is found to approximate the 
full leading-order result quite well as long as the Debye mass is less than the temperature. 



I. INTRODUCTION 



Transport coefficients, such as viscosities, diffusivities, or electric conductivity, character- 
ize the dynamics of long wavelength, low frequency fluctuations in a medium. In a weakly 
coupled quantum field theory, transport coefficients should, in principle, be calculable purely 
theoretically. The values of transport coefficients are of interest in cosmological applications 
such as electroweak baryogenesis [1, 2] and the origin of primordial magnetic fields [3], as 
well as for hydrodynamic models of heavy ion collisions (see, for example, Refs. [4-10], and 
references therein). From a purely theoretical perspective, the evaluation of transport coef- 
ficients also provides an excellent test of our understanding of dynamic processes in thermal 
field theory. 

In a previous paper [11], we performed leading- log calculations of the shear viscosity, elec- 
trical conductivity, and flavor diffusion constants in weakly coupled, high temperature gauge 
theories 1 — that is, neglecting relative corrections suppressed only by powers of the inverse 
logarithm of the gauge coupling, l/log(l/g). Leading-log calculations may be regarded as 
improvements over phenomenological estimates based on relaxation time approximations 
(see, for example, Refs. [14-19]), but no leading- log calculation can be trusted to provide 
even a factor of two estimate in any real application, because the logarithm of the inverse 
gauge coupling, even for electromagnetism, is never all that large. 

At a minimum, one would like to know several terms in the expansion in inverse powers 
of log(l/(?) in order to assess the utility of this asymptotic series. Even better would be a full 
leading- order calculation of transport coefficients, by which we mean an evaluation which 
correctly includes all effects suppressed by powers of l/log(l/g) and only neglects relative 
corrections suppressed by powers of g. This is a feasible goal. Transport coefficients are 
dominantly sensitive 2 to the dynamics of excitations (i.e., quarks and gluons) with typical 
momenta of order of the temperature T. And, as we recently discussed in Ref. [22], it is 
possible to formulate an effective kinetic theory which correctly describes the leading-order 
dynamics of such excitations. The purpose of this paper is to perform a full leading-order 
evaluation of shear viscosity, electrical conductivity, and flavor diffusion constants, and also 
to examine the utility of the asymptotic expansion in inverse powers of log(l/g). We believe 
this work represents the first reliable analysis of transport coefficients beyond a leading-log 
approximation in high temperature gauge theories. 

We begin in Section II by summarizing the relevant definitions, reviewing the form of the 
effective kinetic theory, and sketching how the actual calculation of transport coefficients 
reduces to a linear integral equation which may be accurately solved using a variational ap- 
proach. The evaluation of matrix elements of the linearized collision operator is discussed in 



1 The leading-log result for shear viscosity in pure Yang-Mills theory was first evaluated by Heiselberg 
[12] and closely approximated by earlier work of Baym et al. [13]. A variety of incorrect leading-log 
calculations of transport coefficients are also present in the literature; see Ref. [11] for a discussion of this 
earlier work. 

2 This is true of the transport coefficients under consideration in this paper. It should be noted that bulk 
viscosity requires a separate treatment because it vanishes identically in a scale invariant theory, and so 
requires the correct incorporation of effects which are irrelevant for other transport coefficients. For a 
discussion of bulk viscosity in the case of scalar field theory, see Refs. [20, 21]. 
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Appendices A and B. Our results for the leading-order fermion diffusion constant and shear 
viscosity in QCD are presented in Section III, and the expansion in powers of 1/ log(l/g) is 
discussed in Section IV. We include here next-to-leading-log results for both Abelian and 
non-Abelian theories. We find that a next-to-leading-log approximation reproduces the full 
leading-order results substantially better than one might have expected. Section V presents 
both leading-order, and next-to-leading-log results for the electrical conductivity in QED 
(or in quark plus lepton) plasmas. A short conclusion summarizes our results and comments 
briefly on the relation between our effective kinetic theory treatment and purely diagram- 
matic approaches. Several more technical appendices follow. The first two present details 
of the required 2^2 particle effective scattering rates, and the 1^2 particle effective 
splitting rates which characterize nearly collinear bremsstrahlung and pair production pro- 
cesses in the fluctuating soft thermal gauge field background. A final appendix contains a 
proof that the expansion of leading-order transport coefficients in inverse powers of log(l/<?) 
is only asymptotic, not convergent, and gives a non-rigorous argument suggesting that this 
expansion is not Borel summable. 

II. INGREDIENTS 

Throughout this work, we assume that the gauge coupling (defined at the scale of the 
temperature) is weak, g(T) <C 1. For QCD, this means that the temperature T is asymp- 
totically large compared to Aqcd- We assume that there are no particle masses close to T. 
Zero-temperature masses must either be negligible compared to gT (so thermal self-energies 
are large compared to the mass), or large compared to T (so the particle decouples and may 
be ignored completely). In particular this means we will not consider temperatures just be- 
low thermal phase transitions (or crossovers). We also assume that any chemical potentials 
for conserved numbers are small compared to T. 

A. Definitions 

Transport coefficients characterize a system's response to weak, slowly varying inhomo- 
geneities or external forces. If a weak and spatially uniform electric field E is applied to a 
plasma, the resulting induced electric current is 



where the electric conductivity a is a function of the temperature T, the electric charge e, and 
the particle content of the theory. This constitutive relation is satisfied up to corrections 
involving gradients of the electric field and higher powers of the field strength; strictly 
speaking the conductivity is defined by the above relation in the limit of vanishing frequency 
and wavenumber of an infinitesimal applied electric field. 

In a theory with a conserved global charge (such as baryon or lepton number), the asso- 
ciated charge density n = j° and current density j will satisfy a diffusion equation, 




(2.1) 




(2.2) 
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in the local rest frame of the medium. 3 The coefficient D is called the diffusion constant. 

When one or more diffusing species of excitations carry electric charge, the electric con- 
ductivity is connected to the diffusion constants for various species by an Einstein relation 
(see Ref. [11] for more discussion). If the net number of each species of charge carriers is 
conserved, then 



where the sum runs over the different species or flavors of excitations with e a and D a the 
corresponding electric charge and diffusion constant, respectively. 4 

If the flow velocity of the plasma is not uniform, then the stress tensor (which defines 
the flux of momentum density) will depart from its perfect fluid form. In the local fluid rest 
frame at a point x, the stress tensor, to lowest nontrivial order in the velocity gradient, will 
have the form 



where V is the equilibrium pressure associated with the energy density (T 00 (x)), and the 
coefficients i] and ( are known as the shear and bulk viscosities, respectively. The flow 
velocity u equals the momentum density divided by the enthalpy density (which is the sum 
of energy density and pressure). We will only be concerned with the shear viscosity in this 
paper; the bulk viscosity requires a more complicated analysis, which to date has only been 
performed for a scalar field theory [20, 21]. 

B. Effective kinetic theory 

In a weakly coupled plasma, the energy, momentum, electric charge, and other global 
conserved charges, are predominantly carried by excitations with momenta p ^> gT. Such 
momenta will be referred to as "hard." Weak coupling implies that these excitations are 
long-lived (and hence well-defined) quasiparticles. Thermal corrections to the dispersion 
relations of excitations are order gT in size [i.e., (p ) 2 = p 2 + 0(g 2 T 2 )]. Hence, hard 
excitations are ultrarelativistic. 

The parametrically large separation between quasiparticle lifetimes and the screening 
length of order l/(gT) ensures that the dynamics of these excitations may be reproduced 
by an effective kinetic theory. 5 In other words, one may describe the instantaneous state of 



3 We work throughout in the local rest frame of some arbitrary point x in the system. This frame is defined 
by the Landau-Lifshitz convention of vanishing momentum density (T 0l (x)) = 0. 

4 For a single massless Dirac fermion at vanishing chemical potential, the susceptibility dn a /dfj, a = h T 2 (at 
leading order in the gauge coupling). For quarks, this should be multiplied by 3 (or the number of colors). 
However, the electrical conductivity is dominated by the transport of charged leptons, which diffuse more 
readily than quarks since they only interact electromagnetically. 

5 In particular, the contributions of hard degrees of freedom to the flux of energy, momentum, or other 
global conserved charges are captured by the two-point Wigner functions for the hard degrees of freedom, 
up to sub-leading corrections. As discussed in, for instance, Refs. [20, 21, 23-26], kinetic theory correctly 
describes the time evolution of the two-point function to within 0(g) or smaller corrections. 




(2.3) 



% CP) - V [V, Uj + Vj Ui - | Sij V 1 ui] - C Sij V' u, , 



(2.4) 
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the system by a set of phase space distribution functions f a (x,p,t). The index a will label 
the relevant species of excitations (i.e., gluon, up-quark, up-antiquark, down-quark, down- 
antiquark, etc., for QCD). Formally, the distribution functions represent spatially smeared 
and Wigner transformed non-equilibrium two-point correlation functions. The evolution of 
phase space distribution functions is governed by effective Boltzmann equations of the form 



d d d 

di + Vp 'dx~ + cxt ' dp~ 



f a (p,x,t) = -C a [f], (2.5) 



where v p = dp /dp is the velocity of an excitation with momentum p, F® xt (x,t) is an 
external force acting on excitations of type a (due, for example, to an applied electric 
field), and C a [f] is the collision integral which characterizes the net rate at which species 
a excitations with momentum p are lost (or created) due to scattering processes involving 
other excitations in the plasma. For ultrarelativistic hard excitations, the velocity is a unit 
vector, v p = p, up to irrelevant subleading corrections. 

As discussed at length in Ref. [22] and proposed earlier by Baier, Mueller, Schiff, and 
Son [27] , 6 the collision term needed to reproduce leading-order quasiparticle dynamics 
must contain both 2 <-> 2 particle scattering terms and effective 1 2 particle splitting 
and joining terms. The latter represent nearly-collinear bremsstrahlung and pair produc- 
tion/annihilation processes which take place in the presence of soft (gT scale) thermal fluc- 
tuations in the background gauge field. Hence, the collision terms take the form 

C a [f]=C^[f] + C 1 a ^[f}, (2.6) 

where 

Cr 2 \f](p) = tX-Y. [ \M a c b d (p,k; P ',k')\ 2 (2n)U^(P + K-P'-K') 

^VVa bcd Jkp>k' 

{n P )f b {k) [i±r{p')\ [i±f d m] 

- f c (p') f d (k')[l±r(p)][l±f b (k)}} , (2.7) 
CT 2 [/](P) = 7^J2 rdp'dk'5(\p\-p'-k')^ c (p;p'p,k'p) 



' kp'i 

X 



and 



be 

x 



'r( P ) [i±f\ P 'p)\ [i±nk'p)} - f b ( P 'p)f c (k'p) [i±r(p)]} 
x {f(p) f b (kp)[i±r(p'p)]-n P 'p) [i±np)][i±f b (kp)]y(2.8) 

Here, P = (p°,p), K = (k°,k), etc. denote null four-vectors (so that p° = etc.), v a is 
the number of spin times color states for species a (i.e., 6 for each quark or antiquark and 



6 See also Refs. [28-32] for discussion of the closely related processes which contribute to the leading-order 
photon emission rate. 
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16 for gluons), and J = J d 3 p/[2\p\ (2ir) 3 } denotes Lorentz-invariant momentum integra- 
tion. As usual, upper signs refer to bosons and lower signs to fermions. Ai^(p,k;p',k') 
is the effective two body scattering amplitude for the process ab <-> cd, defined with a rel- 
ativists normalization for single particle states; its square |-M^| is implicitly understood 
to be summed (not averaged) over the spins and colors of all four excitations. Similarly 
7" d (p; p'p, k'p) is the differential rate for an a — > be effective splitting process (or its time 
reverse), integrated over the parametrically small transverse momenta of the participants. 
This rate is likewise understood to be summed over the spins and colors of all three partici- 
pants. The prefactors of l/(4|p|) in the 2 <-> 2 terms are a combination of the l/(2|p|) from 
the relativistic normalization of scattering amplitudes together with a symmetry factor of 
1/2 which corrects for double counting of final or initial states (specifically p', c interchanged 
with k',d). 

The 2^2 matrix elements |-M^| 2 are the usual lowest order vacuum amplitudes, except 
that thermal self-energies must be included on internal gauge boson or fermion lines in t 
and u channel exchange processes [22]. These thermal self-energies are evaluated in the hard 
thermal loop (HTL) approximation, which is adequate for a leading-order analysis. Explicit 
expressions may be found in Ref. [22] and also appear in Appendix A, which discusses the 
evaluation of the 2 <-> 2 collision terms. 

The 1 2 splitting rates encapsulate the effect of near-collinear processes in which 
any number of soft scatterings with other excitations in the plasma occur during the emis- 
sion. The importance of including such multiple soft scatterings is known as the Landau- 
Pomeranchuk-Migdal (LPM) effect. It leads to a two-dimensional linear integral equation 
which must be solved, for each value of \p\ and k, to determine the 1 <-> 2 transition rate 
7^ c . This is discussed further in Ref. [22] and in Appendix B, which describes the evaluation 
of these 1^2 splitting rates. 

C. Linearization 

Given an effective kinetic theory, transport coefficients are calculated by linearizing the 
Boltzmann equation about local equilibrium. The required procedure is discussed in some 
detail in our earlier paper [11]. In brief, one writes each distribution function as a slowly vary- 
ing local equilibrium part }q{x,p) plus a small departure from (local) equilibrium f"(x,p). 
For electric conductivity, /q may equal a homogeneous equilibrium distribution. For dif- 
fusion, /q is a local equilibrium distribution in which the chemical potential coupled to 
the conserved charge of interest varies slowly in space. For shear viscosity, is a local 
equilibrium distribution in which the flow velocity has a non-zero shear. In each case, the 
leading contribution on the left side of the Boltzmann equation (2.5) comes from either the 
convective derivative or external force term acting on /q . Because the collision terms vanish 
identically for any local equilibrium distribution, the leading contribution on the right side 
of the equation is linear in the departure /" from equilibrium. One solves the resulting 
linearized Boltzmann equation for uses the result to evaluate the stress tensor or current 
density of interest, and then reads off the appropriate transport coefficient. 

The convective derivative acting on /q is proportional to the relevant driving term, which 
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we will denote by 

{— Ei , (conductivity), 
Vj Ha ? (diffusion), (2.9) 

(VjWj + VjUi — | 5j 3 -V • u) , (shear viscosity). 

This tensor has I = 1 angular dependence for conductivity and diffusion, and 1 = 1 angular 
dependence for shear viscosity. The departure from equilibrium which solves the linearized 
Boltzmann equation must be proportional to X^..^. It is convenient to express the departure 
in the local rest frame (where fft is isotropic) as 

mp)=p 2 fsm±fs(p)}Xi---jxi. j (p), (2.10) 

which defines xl-j{p)- The resulting linearized Boltzmann equation for the functions xl-j(p) 
may be written compactly in the form 

SIM = (Cxi-iTip) , (2.ii) 

where C is a linearized collision operator defined below. The source term is 

SIM = -Tq a fS( P )[l±fM] h-Ap). (2.12) 

Here q a denotes the relevant conserved charge carried by species a associated with the 
transport coefficient of interest, 

Gem i (conductivity), 

(diffusion), (2.13) 
|p| , (shear viscosity), 

where in the case of diffusion, a is simply a label for the flavor symmetry of interest {e.g. 
quark number or lepton number). 7j...j is the unique t = 1 or £ = 2 rotationally covariant 
tensor depending only on p, 

Pi, i = 1 (conductivity/diffusion), 

The normalization on was chosen so that 

Wp) /«...,-#) = 1 (2.15) 

(which is the reason for the peculiar normalization of in the shear viscosity case). More 
generally, 

I i ... j {p)I i ... j {k)=P e {p-k), (2.16) 

where Pe{x) is the €th Legendre polynomial. Rotational invariance of the collision operator 
(in the local rest frame) implies that 

XlM = I,..j{p) X a {p) (2.17) 
7 



for some scalar function x a (p)- Here and throughout, p = Solving the linearized 

Boltzmann equation means determining x a {p)- 

Let N s be the number of relevant species, and define an inner product on iV s -component 
functions of momenta, 

(/,s)^ 3 £^/^/W(p). (2.18) 

a 

One may show that the linearized collision operator C is symmetric with respect to this 
inner product. It is a positive semi-definite operator, which is strictly positive definite in the 
£=1 and 1=1 channels relevant for diffusion or shear viscosity. Consequently, the linearized 
Boltzmann equation (2.11) is precisely the condition for maximizing the functional 

Q[X] = (x.-r^-j) ~ \ (xi-jXXi-j)- (2-19) 
The explicit forms of the source and collision parts of this quadratic functional are 

d 3 p 



(x,- J .S l ..,)--3 2 Yl J 4&sMP)[^fo(p)]q a X a (p), (2.20) 

(Xi-j, C X i-?) = (Xi-j, C 2 " 2 x,-) + (\,.,. C 1 " 2 \,.., ; ) , (2.21) 

th 

( X ,.,,C 2 ~ 2 x,.,) =^f \M%(p,k;pi,k>)\ 2 (2n)UW(P+K-P'-K') 

V / 5 , , Jpfcp'fc' 



and 



abcd J P k P' k ' 

x / °(p) f b o(k) [l±/ c (p')] [l±tf (*')] 



X 

?3 



'xlM + Xlj(k) - XU{P') - Xt,m] , (2.22) 



x / V) [l±/ >)] [l±/o c (fc)l xV) - x\p) - X c (fc) 



2 



(2.23) 

We have used the normalization condition (2.15) to simplify the expression (2.23). [One 
may also use the identity (2.16) to contract the spatial indices in the 2^2 piece (2.22).] 
The overall coefficients of 1/8 and 1/2 in the collision parts (2.22) and (2.23) are symmetry 
factors which compensate for multiple counting of the same physical process. We have 
dropped the directional information p in the arguments of the collinear splitting functions 
7^ c because, in the linearized theory, these are to be evaluated in the background of the 
isotropic distribution /o(p) and so do not depend on direction. See Appendix B for further 
details. 

Maximizing the functional Q[x] is equivalent to computing an expectation value of the 
inverse of the linearized collision operator, 

Qm&x 2^^"'-?' ^ ^i- j^j ■ (2.24) 
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The value of the extremum directly determines the actual transport coefficients, 



O = | Qr 



D a = iQ 



(2.25a) 

9n - (2.25b) 



1=2, q=\p\ 

The charge susceptibility appearing in expression (2.25b) for D a is 

dn n T 2 



e=l, q=q a \d/J, a 

(2.25c) 



d(j, a 12 



5>aAa(<£) 2 , (2-26) 



where A a is 1 for fermionic and 2 for bosonic species. If n a is the net number density of a 
single fermion flavor, then this susceptibility is | iV c T 2 for a massless Dirac fermion in the 
fundamental representation of SU(iV c ). 

In numerical computations, we will focus in this paper on transport coefficients at zero 
chemical potential, so that the local equilibrium distributions /q (p) appearing in the pieces 
(2.20) and (2.21) of Q[x\ are simply 

/o(p) = ^Vl" (2-27) 

We will primarily discuss the quark number diffusion constant D q and the shear viscosity rj 
in QCD, although we will also report results for the electrical conductivity in a QED plasma 
(with or without quarks). Because we treat all relevant quarks as massless, the diffusion 
constants for individual flavors of quarks in QCD are all the same and equal to D q , which 
is also the diffusion constant for baryon number. 



D. Variational solution 

To maximize the functional Q[x] exactly, one must work in the infinite dimensional space 
of arbitrary functions x{p)- However, as with many other variational problems, one can 
obtain highly accurate approximate results by performing a restricted extremization within 
a well chosen finite dimensional subspace. We will select a finite set of basis functions, 
{(/>( m )(p)}, m — 1 . . . K, and only consider functions x a (p) which are linear combinations of 
these basis functions, 

K 

X a {p) = Y,^n ( t )(m) ^)- (2-28) 

m=l 

(More precisely, we will choose increasingly large sets of basis functions, in order to study 
the convergence of results with the size of the basis.) Restricted to this subspace, the source 
and collision parts of the functional Q[x] become linear and quadratic combinations of the 
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arbitrary coefficients {Xm}i 

(*■■*£■•■;) = J2^ a mS a m , (2.29) 

a,m 

(A,.,- Cx,.,) = EE^^Xn- (2-30) 

a,m b,n 

The matrix elements C^ n may be regarded as forming a square matrix C of dimension 
N S K, while the components {S^} and {Xm} comprise vectors S and x of the same size. 7 
The functional restricted to this subspace is Q[x] = x J S — \x Y Cx- The extremum of 
Q[x] occurs at x = C^S, and 

Qm^ = \x Y S = \S~ T C- l ~S. (2.31) 

Note that the variational estimate Qmax gives a lower bound on the true extremum Q max , 
and a nested sequence of variational estimates must converge monotonically upward to the 
true result. 

The functional Q[x] is most sensitive to the behavior of x a {p) f° r momenta near T. Having 
factored /q [1±/q] out of the first order correction to local equilibrium, as done in Eq. (2.10), 
one may show 8 that the resulting functions x a (p) grow no faster than p e as p — > oo. Also, 
X a (p) vanishes at least as fast as p e_1 at small p. Consequently, one reasonable choice of 
basis functions is 

^ (W) (P) = ^+p/^2 . ™ = 1,-,*- (2-32) 

These functions are not orthogonal; they do not need to be. As noted below, choosing 
strictly positive basis functions improves the accuracy of numerical integrations. These 
basis functions span nested subspaces and, as K — > oo, the basis set becomes complete. 9 



7 As discussed in Ref. [11], when computing transport coefficients in a plasma with vanishing chemical po- 
tentials, charge conjugation symmetry relates and constrains the departures from equilibrium for different 
species. For a QCD-like theory where all matter fields are the same type [e.g., fundamental representation 
fermions), the net effect is that there is only a single independent departure from equilibrium in the C-odd, 
t = 1 channel relevant for flavor diffusion, while there are two independent functions (fermion and gauge 
boson) in the C-even, I = 2 channel relevant for shear viscosity. Hence the actual size of the matrix C 
one needs to deal with is just the basis size K for diffusion, or 2K for shear viscosity. 

8 In the absence of 1 <-> 2 processes, the asymptotic large and small p behaviors are both proportional to p e . 
Because of LPM suppressed 1 <-> 2 processes, the large p behavior is modified to p 2 . However, 2^2 
processes are much more efficient so this behavior only sets in at very large momenta which are irrelevant 
for transport coefficients. The soft region is also modified by 1 <-» 2 processes, leading to p 1 behavior. 
This only occurs for p < mu , and is an accident of our including finite Debye screening masses in 2 <-» 2 
processes but treating 1 <-» 2 processes as perfectly collinear. 

9 For each K, our functions (j)^ are linear combinations of the first K functions in the sequence Xn(p) = 
p e ~ 1 (l+p/T) 1 ~ n , n = 0, 1, . . . . Therefore, as we increase K, the span of our functions strictly increases. 
One may easily show that the functions {xn} form a complete basis in L 2 (R + , dp,) with measure dp = 
w(p)p 2 dp and w(p) any weight function which is real, positive, smooth, bounded, and falls faster than 
p~ 2l ~ 4 as p — > oo. The most natural weight function for our application is w(p) = /o(j?)[l±/o(p)]. 
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that we expect the large K limit of the finite K variational extremum to give the extremal 
value in the full space of allowed \- 

Evaluating matrix elements of the linearized 2^2 collision operator (2.22) requires 
performing an eight dimensional integral (after accounting for energy and momentum con- 
servation). By suitably choosing variables, as discussed in Appendix A, three of these 
integrals represent overall rotations in momentum space and are trivial, and one remaining 
angular integral may be performed analytically. This leaves a non-trivial four dimensional 
integral which must be evaluated numerically for each pair of basis functions. It is computa- 
tionally challenging to perform these numerical integrations both efficiently and accurately. 
Use of a basis of strictly positive, non-orthogonal functions, such as the set (2.32), prevents 
cancellations between different regions of the integration which would otherwise degrade the 
accuracy of the numerically integration. Of course, one does not want to use a basis which 
causes the resulting matrix C to be so ill-conditioned that its inversion in formula (2.31) 
becomes a problem. The choice of the functions (2.32) was motivated by the need to balance 
these two conflicting goals. 

Due to the collinearity of the 1 <-> 2 processes, computing matrix elements of this part of 
the linearized collision operator (2.23) involves only a two dimensional integral. As noted 
earlier, however, evaluating the integrand itself requires the solution of a two-dimensional 
linear integral equation. This evaluation is discussed in Appendix B; one can either convert 
the integral equation to a differential equation/boundary value problem which can be solved 
with a shooting method [32, 33], or one may again use a finite basis variational approach 
[30]. 

E. Thermal masses 

In the next two sections we will present results for SU(iV) or U(l) gauge theories with iVf 
Dirac fermions in the fundamental representation. 10 Because the weak coupling behavior 
of transport coefficients in hot gauge theories is not given by a simple power series in g, 
leading order results cannot be presented just by reporting a coupling-independent value of 
the leading coefficient. Rather, leading order results contain non-trivial dependence on the 
ratio of the effective thermal mass for hard gauge bosons to the temperature, m e fj ig /T, and 
the analogous ratio involving the effective thermal mass for hard fermions, m C Q^/T. These 
masses depend on the gauge coupling and on the particle content of the theory; at leading 
order, 

m e 2 ffig = \ (c A + N { C F Q 5 2 T 2 e|<, (2.33) 

"&,f = iC F g 2 T 2 = 2m 2 F . (2.34) 

As indicated here, the thermal masses for hard excitations, m e g ig and m e ffj , differ by factors 
of \[2 from the more commonly used Debye mass m D or thermal quark mass m F (which 

10 In a theory with chiral (Weyl) fermions, such as the standard electroweak theory (ignoring Yukawa inter- 
actions), each Weyl fermion contributes half as much as a Dirac fermion; so for instance, SU(2) theory 
with 12 left handed doublets behaves the same as SU(2) theory with 6 Dirac doublets. 



11 





40 r 




35 - 




30 - 




25 - 






<* 






20 - 


X 




Q 






15 - 




10 - 




5 - 




- 




FIG. 1: Leading-order value of the quark flavor diffusion constant D q , multiplied by g 4 T, plotted 
as a function of mj)/T. The different curves show the result for SU(3) gauge theory with to 
6 flavors of quarks. The iVf = curve is the result when one artificially neglects scattering of 
quarks off other quarks, and only includes scattering of quarks off gluons (in which case the result 
is independent of Nf ) . 



is the thermal energy of a fermion with zero momentum). We will generally use the latter 
variables below. Here, d-p and g?a denote the dimensions of the fundamental and adjoint 
representations, respectively, while CV and C A are the corresponding quadratic Casimirs. 11 
Note that the ratio of (leading-order) thermal masses, m eS ^/m eS j, is independent of the 
gauge coupling and only depends on the particle content of the theory. 



III. LEADING ORDER DIFFUSION AND SHEAR VISCOSITY IN QCD 

In Figure 1 we plot the fermion flavor diffusion constant D q (multiplied by g A T) as a 
function of mrj/T for QCD with several different values of Nf . Figure 2 shows the shear 
viscosity, multiplied by g A /T 3 , as a function of mv/T for the same set of theories. The results 
in these figures were computed using four basis functions of the form (2.32). Truncation 
error due to the use of a finite basis set is smaller than one part in 10 4 , or smaller than the 
width of the lines in Figs. 1 and 2. Errors dues to the finite numerical integration precision 



11 For U(l), d F = d A = C F = 1 and C A = 0. 

For SU(2), d F = C A = 2, C F = 3/4, and d A = 3. 
For SU(3), d F = C A = 3, C* F = 4/3, and d A = 8. 
For SU(iV), d F = C A = N, C F = {N 2 -l)/{2N), and d A = N 2 -l. 
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FIG. 2: Leading-order value of the shear viscosity rj, multiplied by g 4 /T 3 , as a function of uid/T. 
The different curves show the result for SU(3) gauge theory with to 6 flavors of fermions. 

are also smaller than the widths of the lines. If only two basis functions are used, then errors 
of about 1% result. 

A. Quality of various simplifying approximations 

To illustrate the relative importance of different parts of the collision term, and to examine 
the sensitivity of results to the correct inclusion of thermal self-energies, we plot in Figure 3 
the relative change in the quark diffusion constant of three flavor QCD resulting from various 
(over) simplifications. The different levels of damage follows. 

1. Neglect LPM corrections to the rate of 1 <-> 2 processes. 

2. Neglect collinear 1^2 processes altogether. 

3. Neglect 1 «-> 2 processes and those 2 «-> 2 terms which do not contribute at leading-log 
order. In other words, include only the underlined terms in the 2 <-> 2 matrix elements 
shown in Table III in Appendix A. 

4. Drop all terms which do not contribute at leading-log order, and then replace the 
correct momentum-dependent HTL self-energies in the exchanged propagator by just 
the Debye mass or thermal quark mass. That is, replace 1/t 2 by 1 / (t—m^) 2 in gauge 
boson exchange diagrams, and replace 1/t by t/ (t— m|) 2 in fermion exchange diagrams. 

None of these approximations are correct at leading order in g, but some are far more 
numerically significant than others. As one may see from Fig. 3, neglecting the LPM effect 
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FIG. 3: Sensitivity of the quark diffusion constant in three flavor QCD to various simplifying 
approximations. Each curve shows the ratio of the answer with the indicated approximations 
made, to the full leading order answer. 



makes quite a small change in the results, as does dropping 2^2 terms which do not 
contribute at leading-log order. 12 However, the collinear 1 <-> 2 processes are important at 
the 10% level, while the approximation sometimes used in the literature of replacing the 
bosonic self-energy by just the Debye mass is quite poor and results in errors at the 50% 
level. 



B. Sensitivity to higher-order corrections 

The leading-order results have relative corrections of order g, and hence are only reliable 
when mrj/T is sufficiently small. Of course, for a given level of precision, specifying just 
how small is "small enough" is not possible without knowing more about the actual size of 
sub-leading corrections. 

It should be emphasized that our leading order results do depend on specific choices 



12 Both of these approximations actually overestimate the collision term. The LPM effect suppresses scat- 
terings, so its neglect increases the collinear splitting rate. And, in SU(3) gauge theory, dropping the 
subdominant 2 <-> 2 terms increases the collision rate because the most important neglected term is an 
interference term which is negative in a non-Abelian theory. 
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which were made in defining the linearized collision operator C. We would obtain somewhat 
different numerical results if we had included thermal self-energies on internal lines of dia- 
grams even when this was not strictly necessary, if we had used full one-loop self-energies 
instead of their HTL approximation (valid for momenta small compared to T), if we had 
included thermal corrections to the on-shell dispersion relations in 2 <-> 2 processes, or if we 
had not approximated the nearly collinear 1 <-> 2 processes as strictly collinear. All of these 
effects produce relative changes suppressed by at least one power of g, so we are justified in 
neglecting them. However, handling any of these issues differently can lead to different, but 
equally valid, formally leading-order results. 

Examining the sensitivity of results to changes in the precise definition of the linearized 
collision operator C (which are equally valid at leading-order) is one way to get a "hint" 
as to the likely size of some actual 0(g) sub-leading effects. Fig. 4 shows results for the 
diffusion constant in three flavor SU(3) gauge theory, computed using several different, but 
equally valid at leading-order, definitions of the linearized collision operator. 

We have only considered modifications of the way in which the thermal gauge field self- 
energy is introduced in t or u channel gauge boson exchange diagrams. We focus on these 
modifications because these diagrams numerically dominate the collision term, and because 
such modifications are easy to study. Our "standard" choice is discussed in Appendix A; 
it consists of writing t (or u) channel exchange diagrams as the analogous result for scalar 
quarks, computed with the HTL self-energy included in the gauge boson propagator, plus a 
spin-dependent infrared safe remainder in which the HTL self-energy may be neglected. Our 
first alternative is to instead write t or u channel gauge boson exchanges as the result for 
fermion-from-fermion scatterings, plus a (different) IR safe piece. Expressions for fermion- 
from-fermion scattering (with the HTL gauge boson self-energy included) may be found in 
Appendix B of Ref . [34] . 

The inclusion of hard thermal loop self-energies on t or u channel exchange lines is only 
justified for soft exchange momentum, since otherwise the HTL approximation to the full 
one-loop self-energy is invalid. For large exchange momentum, inserting the HTL self-energy 
is no more correct than neglecting the self-energy altogether, but in this regime self-energy 
corrections are 0(g 2 ) effects. So in a leading-order treatment, one valid approximation is 
to multiply the HTL self-energy by the step function 0(T 2 — Q 2 ), with Q the exchanged 
4-momentum. 13 Alternatively, one may multiply the HTL self-energy by 0(T 2 — q 2 ), with 
q — \q\ the spatial 3- momentum transfer. 

These various possibilities are compared in Fig. 4. The differences between these curves 
are all formally at most 0(g). The width of the band of results, for a particular value 
of itld/T, may be viewed as a guess as to the size of certain types of actual sub-leading 
corrections. The figure also includes the next-to-leading logarithmic approximation from the 
next section, for comparison. One sees that the different, but equally valid, leading-order 
implementations agree to within 15% provided ttid/T < 0.8. This suggests that the leading- 
order results might have about as large a range of utility as one could have reasonably hoped 
for; certainly one should not expect a leading-order weak-coupling analysis to be accurate 
when mo > T. 



We use a ( — h++) spacetime metric, so Q 2 > for a spacelike momentum exchange. 
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FIG. 4: Leading-order results for the flavor diffusion constant D q , multiplied by g T, as a function 
of mn/T in three flavor SU(3) gauge theory. Each curve is computed using a different, but equally 
valid, leading-order definition of the effective kinetic theory. The solid line shows the result of the 
implementation discussed in detail in Appendix A, in which t or u channel matrix elements are 
written as an HTL-corrected scalar quark contribution plus an IR-safe spin-dependent remainder. 
The dotted line is the result of the analogous procedure using fermionic, rather than scalar, scat- 
tering as the template for t channel gauge boson exchange. The dashed lines show what happens if 
the hard thermal loop self-energy is only included in gauge boson exchange lines when the exchange 
momentum Q satisfies Q 2 < T 2 , or alternatively q 2 < T 2 . 

An obvious question is whether there is some natural way to define a unique leading- 
order result which would systematically drop all 0(g) corrections while retaining the full 
dependence on l/lng" 1 . This is directly related to the summability of the asymptotic 
expansion in 1/ In p -1 , and will be discussed in Appendix C. 



IV. EXPANSION IN INVERSE LOGS 



As discussed in Refs. [11-13], part of the linearized collision operator C contains a log- 
arithmic enhancement proportional to ln(T/mD) ~ ln(l/g) arising from t (or u) channel 
2^2 exchange processes with momentum transfer between the Debye screening scale mrj 
and T. One may separate the full leading-order collision operator C into a logarithmically 
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enhanced piece, which we will denote by Cll, and a remainder SC, 



C = C LL + 5C. (4.1) 

Expanding the inverse collision operator appearing in Q max [c./., Eq. (2.24)] in powers of 5C 
will generate an asymptotic expansion of transport coefficients in powers of l/lng -1 . 

To make this separation precise, let C = C/(g 4 T) and note that this rescaled operator 
is dimensionless and depends on g only through the ratios m D /T and m F /T. Equivalently, 
since the (leading-order) ratio mp/mD has a fixed value in a given theory, the g dependence 
of C may be regarded as arising only through dependence on m-o/T. This dependence may 
be isolated by introducing a separation scale satisfying tod C g t < T and splitting the 
relevant t and u channel exchange parts of matrix elements of C into contributions from 
exchange momentum less or greater than q*, 

(xi-j, Cxi-j) = (xi-j, C [q . ,,. Xi-j) + (xi-j, <?[</>?*] Xi-j) • (4-2) 

When q > q*, one may safely set mp, = mp = 0. This contribution is therefore independent 
of mc and mp, up to corrections subleading in g. When q < q*, one may expand in q and 
keep only the first nontrivial order. (This is safe because in all but a g 2 suppressed part of 
the integration domain, all other momenta are large compared to q*.) As far as the q and 
rriD dependence is concerned, the result has the form [11, 12] 



dq f(f) ■ (4.3) 



Rescaling q — > g/mo, the integrand becomes g independent and all remaining g dependence 
is isolated in the upper limit of q^/m^. Simplifying the integrand at q = q*, using the 
assumed scale separation m D <C < T, leads to an explicit form for 

. dC(m D /T) 

A= hm -m D 7^ — ' — -, (4.4) 

m D /T^o am D 

given in Ref. [11]. For a gauge theory with N { Dirac fermions all in the same representation, 
N { d F Cj(3 2 r °° 



(Xi—ji -AXi—j) 



32vr 3 
d A C A (3 3 ml r 



f~dppfi(p)[l+f§(p)] (lx f (p)-x 9 (p)} 2 + [x f {p)-x 9 ip)f) 



(2tt)3^ 
dp Cp / 



o 



dpf 9 M\i+f 9 M\ (p 2 ix 9 (p)? + tm)lx 9 (p)} 2 ) 



^dpc^mi r dpf f m _ f f {p)] U[ X ^ P yf + mi)[ X f 







+ 



P 2 [x f \p)? + t(M)[x f (p)} 2 ) 



(4.5) 



Here x 9 ■, X^ and x^ are the departures from equilibrium for gauge bosons, fermions, and 
anti-fermions, respectively, and primes denote derivatives. The above form assumes that 
departures from equilibrium are fermion flavor independent, as is appropriate for computing 
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shear viscosity or baryon number diffusion. The prefactor of f3 2 m^,/g 2 in two of the terms is 
just a compact way of writing |(Ca + NfCp dp /d&). The first integral comes from fermion 
exchange diagrams whose infrared behavior is regulated by trip, while the latter two integrals 
arise from gauge boson exchange diagrams whose IR behavior is regulated by m D . 

Since the operator A is itself rriD independent, the definition (4.4) implies that the limit 

SC(fi/T) = lim \c(m D /T) - A ln(///m D )l (4.6) 

exists. Therefore the original linearized collision operator, up to 0(g) corrections, may be 
written as 

C = C LL (n) + 5C(n), (4.7) 

where 

ClM = g A TA\n(fi/m D ), (4.8) 

and 

5C(ix) = g A T5C(fi/T). (4.9) 

Here \x is an arbitrary scale which should be chosen to be 0(T) so that 5C(fi) does not 
contain large logs. 

Expanding the inverse collision operator in a geometric series, 

oo 

c- 1 = [C L M + 5C(v)]- 1 ~ (-ITClM' 1 [SCMCvlO*)- 1 ]" , (4.10) 

n=0 

directly gives its asymptotic expansion in powers of [ln(///mD)] _1 ~ 1/ \icig~ 1 . Inserting this 
expansion into Q max = \ (Si-j, C -1 «Sj...j) and using the definitions (4.8) and (4.9) yields the 
inverse log expansion of Q max (and hence of transport coefficients), 

Qmax ~ QM In , (4.11) 

with 

QM ee \ (-1)"- 1 (S t ...„ A- 1 [SC(pi) A' 1 }" 1 ' 1 Si.-.j) . (4.12) 

Except for Qi, the coefficients Q n (^) are /i-dependent. However, the series is formally ji 
independent in exactly the same way that perturbative series in QCD are formally renor- 
malization point independent even though individual n-loop contributions do depend on the 
renormalization point. 

To avoid the presence of large logarithms in the coefficients Q n (n), one must choose // to 
be 0(T) parametrically, but the exact coefficient is not uniquely prescribed. One somewhat 
natural choice is to select the value /i* for which (/-**) vanishes. This can be termed the 
fastest apparent convergence (FAC) choice of \x at next-to-leading log order. Given results 
for Qi and (^(/-O at some other value of //, 

//* = // exp[-Q 2 (n) /Qi] . (4.13) 
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flavor diffusion 


shear viscosity 


group 




Di 




m 


»*/T 


U(l) 


1 


47.089 


2.469 


188.38 


5.007 


2 


30.985 


3.013 


120.28 


4.418 


SU(2) 


6 


21.283 


3.123 


200.533 


2.927 







16.060 


2.699 


27.126 


2.765 




2 


12.999 


2.887 


86.47 


2.954 




3 


11.869 


2.949 


106.66 


2.957 


SU(3) 


4 


10.920 


2.997 


122.96 


2.954 




5 


10.111 


3.035 


136.38 


2.947 




6 


9.414 


3.065 


147.63 


2.940 




oo 


136.76/7V f 


3.155 


274.83 


2.733 



TABLE I: Values of the leading- log coefficient D\ and rj\ together with the value of fJ,*/T, for the 
case of fermion flavor diffusion and shear viscosity in theories with the indicated gauge group and 
Nf Dirac fermions in the fundamental representation. The large Nf result in the last line is from 
Ref. [34], and shows that Nf = 6 is still a long ways from the large Nf limit. The values shown for 
Nf =0 flavor diffusion represent the results one would obtain if diffusing quarks could only scatter 
off gluons, and not off other quarks. 



To evaluate the coefficients Q n {y) of the expansion (4.11), we use exactly the same finite 
basis set approach described in section II D. The linear operators A and 5C are replaced by 
their matrix representations A and 5C in the finite basis set, and 

4nax ~ -t= QM Un , (4.14) 

with 

QM = |(-1)" _1 S T A- 1 ^C(fi) i- 1 ] S . (4.15) 

The limit (4.6) defining 5C(fi) is performed numerically by evaluating each matrix element 
of C at several small values of mj), subtracting off the leading log piece (whose matrix 
elements are easy to evaluate), and then extrapolating to vanishing Debye mass. We find 
that this extrapolation is quite well behaved, although the numerical integrals at small m D 
or mp become rather demanding. 

Using the basis functions (2.32), we previously found [11] that the fractional difference 
between Q\ and its finite basis approximation Q\ is less than 10 -5 with 4 basis elements 
and less than 1CT 6 using six. The higher Q n are more sensitive both to basis size and to 
numerical integration errors: Q2 can be reliably determined with 4 basis functions, and Q3 
and Q4 can be found with reasonably small errors using 6 to 8 basis functions, but higher 
moments become rapidly more difficult to evaluate, showing poor convergence with basis set 
size and high numerical integration error sensitivity when the basis sets become very large. 
Consequently we have been unable to go very deep in the Q n series. 

Table I shows results for //* and the first coefficient of the inverse log expansion, for the 
case of fermion flavor diffusion and shear viscosity in theories with various gauge groups 
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and the indicated number JVf of Dirac fermion flavors in the fundamental representation. 
Specifically, we show the first coefficients D x and 771 of the series 14 



1 00 / 



n=l 



JJ_ 

m D 



and 



^ ~ Vniv) fin— ) 



(4.16) 



(4.17) 



Together with the corresponding values of /i*/T, these numbers determine the next-to- 
leading-log (NLL) approximation to the respective transport coefficients, 



-Dntj, — 



9 4 T 



D 1 



and 



T 3 
T 



5 



ln(/i*/m D ) 
ln(/i*/m D ) 



(4.18) 



(4.19) 



Note that in SU(3) gauge theory, /i*/T is quite close to 3 for both transport coefficients, 
regardless of the number of fermion flavors. 

We have computed further terms in the inverse log expansion in the case of SU(3) gauge 
theory with 3 fermion flavors. For flavor diffusion, we find 



DM = 2.436(2) , DM = -0.11(1) , DM = 1.7(1) , 
while for shear viscosity 



V M = 27(1) , 



77 4 (//*) = 6(5) , 



r)M = 100(100) . 



(4.20) 



(4.21) 



[And D2 (/•**) — = ?72 (/-**)> ^ our definition of //*.] The third-order coefficients D$(n*) and 
773 (/J*) are roughly one quarter the size of D\ and 771, respectively. The next order coefficients 
D^/j*) and ^(/i*) are yet smaller, but subsequent coefficients appear to grow. 

The behavior of the expansion in inverse powers of ln(/j*/mD), truncated at second or 
third order, is compared to the full leading-order result in Figure 5, for the case of flavor 
diffusion in three flavor QCD. As the figure makes clear, the next-to-leading-log result is 
remarkably close to the full leading order result out to m D /T = 1, but going beyond second 
order in the inverse log expansion has very little practical utility. 

Appendix C discusses the asymptotics of the inverse log expansion and proves that this 
expansion is, in fact, a typical asymptotic expansion with zero radius of convergence. An 
argument is also given suggesting that this asymptotic expansion is not Borel summable. As 
is well known, the presence of singularities in the Borel transform on the positive real axis 



The diffusion constant D and shear viscosity 77 are related to their respective Q ma x's as shown in Eq. (2.25). 
Explicitly, 77 = ^ Qmax and D = |<5max/ (Nf T 2 ) with Q^ D \ in this normalization, being associated with 
total quark number. Hence, the inverse log expansion coefficients in Eqs. (4.16) and (4.17) are related to 
the previous coefficients Q n {p) via D n (y) = (p)/(N{ T 2 ) and rj n (^i) = j^Q^\/j,)/T 4 . 
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FIG. 5: Leading-order results for the flavor diffusion constant (multiplied by g^T) in three flavor 
QCD compared to the expansion in inverse powers of ln(/i»/mD) truncated at second and third 
order. Right panel: zoom-in on the small m£)/T region, showing that the third order truncation 
of the inverse- log expansion can be an improvement over the second order result, but only for 
m D /T < 0.2. 

generates ambiguities in the inverse Borel transform. In the typical case of a power series 
in g 2 , this irreducible ambiguity is non-perturbative, behaving as exp(— c/g 2 ) where c is the 
location of the singularity of the Borel transform nearest to the origin on the positive real 
axis. (See, for example, Refs. [36, 37].) In the present case where the expansion parameter 
is an inverse log of the coupling, a singularity in the Borel transform instead indicates an 
inherent ambiguity in the Borel sum of the asymptotic series which is a power of coupling. 
The estimate of appendix C suggests that this is an 0(g 5 ) ambiguity. 



-Leading-order 
NLL FAC 

- NLL FAC + NNLL term - 



V. ELECTRICAL CONDUCTIVITY 



We have also computed the electrical conductivity for a high temperature plasma of 
leptons, or leptons plus quarks. As mentioned earlier, the electrical conductivity is related 
to the diffusion constants of charged species via the Einstein relation (2.3). In a plasma 
of leptons plus quarks, we make the same a| M -C a 2 approximation used in previous work 
[11, 35]: we neglect the electric current directly carried by quarks, and only compute the 
charged lepton diffusion constant. Because quarks undergo efficient QCD scattering (as 
compared to QED scatterings) , their departure from equilibrium in the presence of an electric 
field is negligible compared to that of charged leptons, and hence so is their contribution 
to the electric current. This approximation amounts to the neglect of relative 0(a| M / 'a 2 ) 
corrections to the conductivity. Quarks remain relevant, however, as excitations off of which 
charged leptons can scatter. Of course, for plasmas containing quarks we still require a s <l, 
so that quarks may be treated as nearly free massless excitations. 
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leptons 


quarks 


Eg 2 


0"i 


M*/T 


a x e 2 /T 
NLL FAC leading order 


e 





1 


15.696 


2.470 


5.928 


5.9700 


e,fi 




2 


20.657 


3.013 


8.262 


8.2996 


e,n 


u, d, s 


4 


12.287 


3.268 


5.498 


5.4962 


e,n,r 


u, d, s, c 


19/3 


12.520 


3.306 


6.208 


6.1756 


e,n,r 


u, d, s, c, b 


20/3 


11.972 


3.306 


6.013 


5.9769 



TABLE II: Electrical conductivity in plasmas containing the indicated types of leptons and quarks. 
Each entry is relevant for temperatures such that the listed species are much lighter than T while 
all other leptons or quarks are much heavier than T. Relative corrections of order a^ M /a 2 are 
neglected; see text. The last two columns compare the next-to-leading-log (NLL FAC) approx- 
imation with the full leading-order result, both evaluated at the physical value of tod (using 
"em = 1/137.04). Clearly, the NLL FAC approximation works very well for QED. 



For simplicity we have only analyzed physically relevant combinations of leptons and 
quarks. We have evaluated the electrical conductivity o at next-to-leading-log order (NLL 
FAC) as well as at full leading order. In the next-to-leading-log form, 



T 

C"NLL — —z 



_ln(//*/m D ). 

the Debye mass now refers to the inverse QED screening length given by 



ml 



±e 2 T 2 



(e,<) 



(5.1) 



(5.2) 



qi is the charge assignment of a given species, and the sum runs over all Dirac fermions. 
Instead of presenting plots showing the conductivity as a function of oem, we have simply 
set a E M = 1/137.04. 15 Our results are presented in Table II. For this quite small value of 
coupling, one sees that the NLL FAC treatment agrees with the full leading order result to 
better than 1%. 



VI. CONCLUSION 



We have performed complete leading-order calculations of shear viscosity, electrical con- 
ductivity, and fermion diffusivity in QCD and QED. "Leading-order" means that all ne- 
glected effects are suppressed by one or more powers of the gauge coupling g(T). To our 
knowledge, this is the first time any transport coefficient has been evaluated with leading- 
order accuracy in a high temperature gauge theory. Due to the presence of Coulomb log- 
arithms arising from small angle scattering, the coefficient of the leading power of g(T) is 
not a simple number, as in scalar theories, but rather is a non-trivial function of ln^" 1 ). 



We use this many significant digits, ignoring the running of the coupling, for no reason other than to show 
the precision of the numerics and to compare different results. 
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Hard off-shell 

Soft, spacelike, gauge boson, HTL resummed 
Hard on-shell, resummed with diagrams of form 




etc. 



FIG. 6: Typical diagram needed in the leading-order evaluation of the shear viscosity in QCD. The 
crosses at the left and right denote (stress tensor) insertions. 

Leading-order results for transport coefficients may themselves be expanded in powers 
of l/ln((yf _1 ). We have explicitly computed both leading and first sub-leading terms for 
shear viscosity and quark diffusivity in U(l), SU(2), and SU(3) gauge theories with various 
numbers of fermion fields (as well as several more terms for three flavor QCD). For QCD, 
the next-to-leading log result (with the sub-leading term absorbed by suitably shifting the 
scale inside the leading log) was found to be remarkably close to the full leading-order 
result as long as ttid/T < 1. This is a much larger domain of utility than one might have 
expected. For these transport coefficients, we also find that only roughly 10% errors are 
made if one neglects near-collinear 1 «-> 2 particle splitting processes, which are considerably 
more difficult to analyze than 2^2 particle scattering processes. (However, it should be 
noted that some transport coefficients which we have not analyzed, such as bulk viscosity, 
depend primarily on particle number-changing processes and so may be expected to depend 
essentially on I <-> 2 processes.) 

Because the expansion in inverse powers of m(g _1 ) is only asymptotic, not convergent, as 
demonstrated in Appendix C, we are not able to give a unique, unambiguous prescription for 
separating leading-order contributions from higher-order effects. As discussed in Appendix 
C, it appears that the inverse log expansion is not Borel summable, which would imply that 
no clean separation is possible. In practice, this means that any specific calculation yielding 
results of leading-order accuracy will necessarily include some contributions from higher- 
order effects. However, our examination of several different prescriptions for computing 
leading-order results suggests that this is not a significant issue for mo ^ 0.8 T. 

Our tool for studying transport coefficients has been kinetic theory, specifically the effec- 
tive kinetic theory presented in our previous paper [22]. As originally shown by Jeon [20], 
in the context of weakly-coupled relativistic scalar theories, it is also possible to compute 
transport coefficients diagrammatically starting from the appropriate Kubo formulae involv- 
ing current- current or stress-stress correlators. Such a diagrammatic approach amounts to a 
complicated way to derive the appropriate linearized Boltzmann equation specialized to the 
particular symmetry channel of interest. For gauge theories, this diagrammatic approach 
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has been applied, only at leading logarithmic order, to the electrical conductivity by Valle 
Basagoiti [38]. Trying to use a diagrammatic approach for a complete leading order calcula- 
tion would be an enormously more difficult task. For instance, a typical diagram which we 
believe contributes at leading order to the (T^Tjy) correlator, needed for the shear viscosity, 
is depicted in Figure 6. Note that the hard (nearly) on-shell propagators require self-energy 
resummations, illustrated at the bottom of the figure, which go far beyond the HTL approx- 
imation. The one and two loop self-energy contributions shown account for scattering via 
2^2 processes, while the very complicated self-energy diagram represents one contribu- 
tion to the effective 1 2 splitting process; the reason it needs so many loops is that this 
process can involve any number of soft scatterings off of other particles in the plasma [the 
number of such scatterings is summed over by Eq. (B3)]. The complicated "cross-rungs" in 
the upper diagram are the result of opening up one of the lines in any one of the self-energy 
contributions. For more discussion of this point, see Ref. [39]. 

An interesting problem for the future is to understand the accuracy of leading-order 
calculations of transport coefficients by calculating higher order effects explicitly. In the 
case of the QCD free energy it is known that, beyond the ideal gas result, the perturbative 
expansion in powers of g(T) is quite poorly behaved [40-42], except for unrealistically small 
values of g(T). (Specifically, T must substantially exceed the Planck scale.) Does this same 
unpleasant behavior apply to transport coefficients? At the moment, the only known test 
case is a many flavor limit of QCD [34] , where the leading-order result (as well as the next- 
to-leading-log approximation thereto) is quite successful — its accuracy is comparable to the 
renormalization point sensitivity. It would be useful to know if this holds more generally. 
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APPENDIX A: 2 2 MATRIX ELEMENTS 

The matrix elements for all 2 <-> 2 particle processes in a QCD-like theory, neglecting 
thermal self-energy corrections, are listed in Table III. These matrix elements arise from 
the diagrams shown in Fig. 7. Terms in Table III with underlined denominators are suffi- 
ciently infrared sensitive that thermal self-energy corrections must be included, as discussed 
in Ref. [22] . Singly-underlined denominators indicate IR sensitive contributions arising from 
soft gauge boson exchange, while double-underlined denominators indicate IR sensitive con- 
tributions from a soft exchanged fermion. 
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(J) 



FIG. 7: Diagrams for 2 <-> 2 processes needed at leading order in the coupling. Leading-log 
calculations require only the squares of diagrams A—E. Next-to-leading-log, or full leading order 
calculations require evaluating the squares of diagrams A— E with HTL self-energies inserted on the 
internal lines, and then adding the (undressed) squares of diagrams F—J, as well as the interference 
terms between diagrams with the same initial and final states. 



ab <-> cd 



MfX/g 4 



qm <-»■ <?i<?2 , 

qm <-»• gi92 , 
<M2 <-»■ M2 

gigi <->• gift , 
Mi <-»■ gift 

gift <-»■ gift 
gift <-»■ g2g2 

gift ^55 

gi 5 <-> gi 5 , 
ft 5 <-> ft 5 

55 55 



4 / / £ 2 +« 2 



4C| /V + u 2 . s 2 + t 2 



d A 



t 2 



+ 



8 



4 Cg /V + u 2 . t 2 + u 2 
dA 



+ 



+ 16d F C F |^C F - 
+ 16d F C F ( C F - 





S 2 










2 , 


I- 



d 2 C 2 /t 2 + u 2 




8 dp C F Ca 



-8d F Cl \ - + - ) + 8d F C F C A 



16 d A C\ ( 3 - —y x ~ 



t 2 + u 2 



s 2 + u 2 



t 2 



TABLE III: Squares of vacuum matrix elements for 2 <-> 2 particle processes in QCD-like theories, 
summed over spins and colors of all four particles. q± and qi represent fermions of distinct flavors, 
qi and q~2 are the associated antifermions, and g represents a gauge boson. Note that the process 
gig2 Q1Q2, for example, appears 2N{ (iVf —1) times in the sum ^2 abcd over species in the linearized 
collision operator (2.19), while qiqi <-> qiqi and gift <-> 55 each appear 47Vf times, (75 <-> (75 appears 
just once, etc. 
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1. Self-energy corrections to matrix elements 



Fermion self- energy 



When 2 <-> 2 particle processes involving t (or -u) channel fermion exchange are computed 
using free propagators, the resulting squared matrix elements (shown with double-underlined 
denominators in Table III) generate logarithmic infrared divergences in the collision term 
(2.22). This logarithmic infrared sensitivity is cut off by the inclusion of the retarded ther- 
mal self-energy S(Q), so that the internal fermion propagator appearing inside the matrix 
element M. is [Q — ^(Q)] -1 . Its conjugate, the advanced self energy S*(Q), appears in 
M*. Since the exchange four-momentum = (u,q) is spacelike, and the thermal self- 
energy is only relevant when Q is soft, we only need the self-energy for spacelike momenta, 
q = \q\ > in the hard thermal loop limit, 5 < T. In this regime the self-energy was 
originally evaluated by Klimov [43] (and independently by Weldon [44]), and is given by 



£°(Q) 
HQ) 



~2q~ 



In 



q + u> 
q — uj 



— 171 



mi 



In 



q + uj 



UJ 



ITX 



(Al) 
(A2) 



Here m F = a/Cf/8 gT is the (leading-order) "fermion thermal mass," equal to the thermal 
energy of a fermion at zero momentum. 

For the process fg — > fg, the net effect of the inclusion of the fermion self-energy is to 
make the replacement 



-|tr[f ^Qrflv&lA _ -4 Re[(P ■ Q)(P' ■ Q*)} + t Q ■ Q* 



u 



\Q-Q\ 



\Q-Q\ 



(A3) 



[with (— +++) metric convention] in the contribution shown in Table III. Here P M = (\p\,p) 
and P /fl = (\p'\,p') are the incoming and outgoing fermion 4-momenta, respectively, and 
= P 1 * — K'^ — Yy-{P—K') with K' the outgoing gauge boson 4-momentum. 
The analogous replacements needed in the // <-> gg squared matrix elements are 



and 



u 
1 

t 
u 



4Re[(P ■ Q)(K ■ Q*)] + sQ - Q* 



\Q-Q\ 

4Re[(P ■ Q)(K ■ Q*)\ + sQ - Q* 



\Q-Q\ 



QV=pV-pili--JlH(P-pi) 



Qv=Pn-K'»-W(P-K') 



(A4) 



(A5) 



where P and K are the incoming fermion momenta, and P' and K' the outgoing gauge 
boson momenta. 



Gauge boson self- energy 

Processes involving t or u channel gauge boson exchange require inclusion of the thermal 
gauge boson self-energy on the internal propagator to cut off the infrared sensitivity of these 
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processes. Because the self-energy only matters when the exchange momentum is soft, one 
may exploit the fact that soft gluon exchange between hard particles is spin-independent 
(to leading order) [31]. If one separates the IR sensitive matrix elements (those with singly- 
underlined denominators in Table III) into (i) the result one would have with fictitious scalar 
quarks plus (ii) a spin-dependent remainder, then all the IR sensitivity resides in the first 
spin-independent piece. This is the only piece which must be recomputed with the thermal 
self-energy included. For the t-channel exchange terms, this amounts to using the exact 
identities 



s 2 +u 2 



2 ~ 2 



i (s-uf 



t 2 



su 



4 f 2 



(A6) 



and then replacing 



[s-u 



t 2 



\D{P-P')^(P+PY(K+K') 



l\v\2 



(AT) 



where D(Q)^ U is the retarded thermal equilibrium gauge field propagator, evaluated in the 
HTL approximation. 

The HTL result for the above replacement (A7) does not depend on gauge choice. One 
convenient choice is Coulomb gauge, where [45] 



-1 



D 00 (uj,q) = 

Dij(uj,q) = 
D 0i (uj,q) = D i0 (u,q 



q 2 + U 00 (u,q) ' 
QiQj 



q 2 - uj 2 + n T (co>, q) ' 



0. 



The equilibrium transverse and longitudinal gauge boson self-energies are [46, 47] 

,2 , , („2 , ,2\ 



U T (u,q) 
Tloo(u,q) 



m 



D 



ml 



lu* u(q 2 —uj 2 ) 

v + V 



In 



q + oj 
q — uj 



l 7T 



1-^ 

2q 



In 



q + uj 



— ITT 



,1 ~ U , 

where we have assumed < q, which is the only case of relevance. 



(A8) 

(A9) 
(A10) 

(All) 
(A12) 



2. Integration variables 

Since all external particles are to be treated as massless, the domain of the phase space 
integrations appearing in matrix elements of the linearized collision operator (2.22) are 
the same for all 2 <-* 2 processes. One rather straightforward method for doing the multi- 
dimensional numerical integration is to use an adaptive Monte Carlo integrator. This can 
give reasonable accuracy at a tolerable investment of computational effort, but for the highest 
accuracy it is preferable to use nested one-dimensional adaptive Gaussian integration. 

In order to handle efficiently the infrared-sensitive terms in t and u channel processes 
which give rise to leading-log contributions, it is useful to pick the exchange momentum q 
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and energy uj as two of the integration variables. When doing nested adaptive quadrature 
integrations, it is especially advantageous to choose integration variables in a manner which 
allows one to perform analytically as many of the integrations as possible. In particular, 
it is convenient to use different parameterizations for terms in Table III having different 
denominators. For terms having denominators of t (or t 2 ), the t channel parameterization 
described below allows all but four integrations to be done analytically. And similarly, the u 
and s channel parameterizations described below allow all but four integrations to be done 
analytically for terms with denominators of u (or u 2 ), or s (or s 2 ), respectively. The constant 
term in the gg <-> gg matrix element can be handled using any of these parameterizations. 
The only other terms in Table III are those involving s 2 /tu or u 2 / st. These can be reduced 
to the previous cases by rewriting 

s 2 s s u 2 u u 

tu t M ' St S t 

In what follows, our convention for labeling momenta in 2 2 processes is that P,K <-> 
P 1 , K' . 



(A13) 



t channel parameterization 

For terms containing t = — (P' — P) 2 in the denominator, it is convenient to use the spatial 
delta function in (2.22) to perform the k' integration, and to shift the p' integration into an 
integration over p'—p = q. The angular integrals may be written in spherical coordinates 
defined such that the z axis is in the direction of q while p lies in the xz plane. This yields 

^3 roc nl p2w j 

(Xi-j,C 2 " 2 Xi~j) = TT^eYl / q 2 dqp 2 dpk 2 dk / dcos9 pq dcos9 kq / d<f> 
V J J (4tt) 6 ^Jo 7-i Jo pkp'k' 



x \M? d \ 2 5(p+k-p'-k') fS(p) f b (k) [1±/ V)1 \l±f^k')] 

x [xl,(p) + Xljik) - xl,(p') - xt, ■(*')] 2 , (A14) 

where p, k, and q denote to the magnitudes of the corresponding three-momenta (not the 
associated 4-momenta), p' = \q + p\ and k' = \k — q\ are the magnitudes of the outgoing 
momenta, is the azimuthal angle of k (and k') [i.e., the angle between the p,q plane and 
the k,q plane], and 9 pq is the angle between p and q (so cos8 pq = p ■ q), etc. 

Following Baym et al. [13], it is convenient to introduce a dummy integration variable u, 
defined to equal the energy transfer p' — p, so that 

/oo 
du 5(u + p - p') 8{u - k + k') . (A15) 
-oo 

Evaluating p' — \p + q\ in terms of p, q, and cos9 pq , and defining t = u 2 — q 2 (which is the 
usual Mandelstam variable), one finds 

5(uj + p-p') = P-5(cos9 pq --- 7 ^-)e(u;+p), (A16) 
pq \ q 2pqJ 

S(u-k + kT) = ^( cos ^-^ + ^) (^-^)' ( A17 ) 
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with G the unit step function. The cos 6 integrals may now be trivially performed and yield 
one provided p > |(g — u), k > ^(q + cu), and \u\ < q; otherwise the delta functions cannot 
both be satisfied. The remaining integrals are 

[Xi-j,C 2 " 2 Xi-j) = jt^eY] / dq du dp dk dxj> 
V } W b ^ d Jo J- g Jz? Jzp Jo 

x \M? d \ 2 r {p) m [i ± /tip')] [i ± f*(k')] 

x [xl.M + Xljik) - Xl-M) - xUik')] 2 , (A18) 

with p' = p + uj and k! = k — cu. For evaluating the final factor of (A18), we use the 
relationship (2.16) that 

h...,(p) h.^k) = P e (cos6 pk ) . (A19) 

One therefore needs expressions for the angles between all species, as well as the remaining 
Mandelstam variables s and u, which may appear in A4 2 . They are 

s = ^{[(p + p')(k + k') +q 2 ] - cos ^(4pp' + t) (4kk' + t) } , (A20a) 
u = -t-s, (A20b) 

and 



COS 6 pq 


UJ t 

q 2pq ' 


cos 9pi q 


UJ t 

q 2p'q ' 


(A21a) 


COS 9 k q 


UJ t 

q 2kq ' 
1 + 2pp' ' 


COS0 k >g 


UJ t 

" q + 2k 'q ' 
+ 2kk' ' 


(A21b) 


cos 6p p i 


cos 9 kk , 


(A21c) 


cos 6 p k' 


1 U 
+ 2pk' ' 


cos 6 p i k 


1 U 
+ 2p'k ' 


(A21d) 


cos 6 pk 


1 2pk ' 


cos 6 p i k i 


= 1 5 . 
2p'k' 


(A21e) 



For terms in which only t appears in the denominator of the matrix element, the inte- 
gration can be easily performed analytically, leaving four integrals which must be evaluated 
numerically. 

u channel parameterization 

For terms in which u = —{K 1 — P) 2 appears in the denominator, exchanging p' and k' 
in the t channel parameterization provides the natural choice of variables. 

s channel parameterization 

For terms in which s = — (P + K) 2 appears in the denominator, one may use the spatial 
delta function in Eq. (2.22) to perform the k' integration, and then shift the k integration to 



29 



an integral over q = p+k, the total incoming spatial momentum (and the momentum on the 
internal propagator in s channel exchange processes). Again choosing spherical coordinates 
so that q lies on the z axis and p lies in the xz plane, the 2^2 contribution (2.22) becomes 

[Xi-jiC 2 ~ 2 Xi-j) = J^eYI J Q 2 dqp 2 dpp' 2 dp' J d cos 9 pq d cos 9 p > q J d( P pkp , k , 

x \Mf d \ 2 5(p+k-p>-k>) fS(p) f b (k) [1 ± r {p>)} [1 ± f^k')] 

x [xl.jip) + Xlj(k) - Xl.jip') ~ Xt.jm} 2 . (A22) 

where now k = \q — p\, k' = \q — p'\, and (f> is the azimuthal angle of k (and k'). Introducing 
the total energy u via 

poo 

5(p + k-p- k!) = duo5{uo-p-k) 5{uo - p - k!) , (A23) 
Jo 

and defining s = uj 2 — q 2 (which is the usual Mandelstam variable), one finds 

5{uj-p-k) = — s(cos9 pq - - + J-}e(u-p), (A24) 
pq \ q 2pq J 

S( u -j/-kf) = ^s(cose p/q -- + -^-)e(u;-p > ). (A25) 
pq \ q 2p'q J 

Integration over cos9 pq and cos9 p i q yields unity provided q < u, \2p—u\ < q, and \2p'—u\ < q 
(and zero otherwise). Therefore, 



abed " u " u " " — 

x \M? d \ 2 fS( P )f b (k) [1±JW)] \l±f d {k')] 

[Xl,(p) + Xl,(k) - xl,(P') ~ Xtjik')} 2 , (A26) 



x 



with k = uj — p and k' — uj — p' . The other Mandelstam variables are 

t = | [{p - k)(p - k') - q 2 ] + cos <f>y/(Apk - s)(Ap'k' - s)} , (A27a) 
u = -s-t, (A27b) 
and the angles between q and the external momenta are 

cos6 pq = - - , cos 0p> q = — - , (A28a) 

q 2pq q 2p'q 

Eqs. (A21c)-(A21e) for the angles between external momenta still hold. For terms in which 
only s appears in a denominator, the integration can be easily performed analytically, 
leaving four integrals to do numerically. 
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APPENDIX B: MATRIX ELEMENTS FOR 1 2 PROCESSES 



In this appendix, we will review the integral equations that determine the splitting/joining 
functions 7^, for effective 1 <-> 2 processes, appearing in our various formulations (2.8) and 
(2.23) of the Boltzmann collision term. These equations are summarized in Ref. [22], but 
here we will make a number of simplifications applicable to the problem at hand. 

The only 1 2 processes with 0(g 4 T 4 ) total rates per unit volume involve hard, collinear 
creation or destruction of a gauge boson. In QCD, the relevant processes are q <-> qg, 
q <-> qg, g <-> qq, and g <-> gg. By "hard" we mean that the outgoing states each have O(T) 
energy, and by "collinear" we mean that the angles between the three external momentum 
vectors are all 0(g). [These parametrically small opening angles are ignored when evaluating 
distribution functions in the 1 <-> 2 piece of our collision term (2.8), which only makes a 
sub-leading error in the evaluation of transport coefficients.] 

The duration of such processes (also known as the formation time of the gauge boson) 
is parametrically 0(l/g 2 T), which is the same as the mean free path for hard particles 
to undergo soft scattering (with momentum exchange of order gT). For this reason, one 
must sum all possible number of soft scatterings with other excitations during the emission 
process. (In realistic theories, at least one such soft scattering is required to allow for energy- 
momentum conservation in the 1^2 process.) This summation can be implemented by 
an integral equation which must be solved numerically. For a complete discussion, and a 
derivation of the initial integral equation presented below, see Refs. [22, 29, 31]. 



1. The integral equation for 7^, 

For non-Abelian gauge theories such as QCD, the splitting/joining functions 7ft c (p';p, k) 
for particle types a <-> be with momenta p' <-> pk are given in equilibrium by [22] 

/2 2 

7 9 W; p. k ) = 7t>'; p, k) = p -j^% Up', p, k) , (Bia) 

P p K 



l 9 qq(p';P,k) = f 2 2 /3 ?<l(k,-P,P') , (Bib) 

At p A p 6 



where 



p/4 _|_ pi _|_ £4 

l 9 gg (p';P,k) = pl3p3k3 F g (p',p,k), (Blc) 

*'W>P>Q = %§$ J-^2h.ReF s (h;p',p,k) (B2) 

and a = ^ 2 /(47r). The function F s (h;p',p, k), for fixed given values of p', p, and k, depends 
on a two-dimensional vector h which is related to transverse momentum during the splitting 
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process. F s is the solution to the linear integral equation 

2h = i5E{h;p',p, k) F s (h;p' ,p, k) 
d 2 q A 



+9 2 J ^ A ^) {(<?• ~ \Ck) [F s (h;p',p,k) - F s (h-kq ± ;p',p,k)} 



+±C A [F s (h;p',p,k) - F^h+p'q^-p'^k)] 
+ ±C A [F s (h;p',p,k)-F s (h-pq ± ;p',p,k)} \ , 



where 



and 



5E(h;p',p,k) 



A{q L ) 



m 



cff,g 



2k 



+ 



m, 



eS,s 



_ _ m eS,s h 2 

2p 2p 2pkp' 



d £{A'{Q)[A-{Q)r)) 



(B3) 



(B4) 



(B5) 



Here Q = (q°,q ± ,q z ), A~ = A - A z , and (A~ (Q)[A~ (Q)]*} is the thermal Wightman 
correlator evaluated in the hard-thermal-loop approximation. Explicit formulas for this 
correlator in equilibrium may be found, for example, in Ref. [31]. However, we will make use 
here of the wonderful simplification found by Aurenche et al. [33] showing that the integral 
(B5) has a remarkably simple form, 



A(q J 



T 



+ m 



D 



(B6) 



2. Variational solution 

One way to solve the integral equation for F s is to use a variational method, similar to the 
method used in the main text for the Boltzmann equation and implemented in the Abelian 
case for 1 <-> 2 processes in Ref. [30]. Some further simplifications are then possible, as we 
shall describe below. The variational formulation is 

J~ i sip' ' 1V1 k) = ^7r) 3 (^ s ) cxtrcmum ' (B7) 

with 

Q s [F] = Re[(2h,F^) -^F,(t5E + IC s )F^ , (B8) 

(Q S )extremum = ^Re^h, (i5E + K 8 ) ~W) , (B9) 

where in this context 

(f,9) = J ^f( h )-9(h). (BIO) 
The analog JC S of the collision operator is given by 

(F,/C S F) =<7 2 {(C.-±Ca) (F,K(k)F) +\C A [F,K{-p')F) + \C a (f ,K(p) f)}, 

(BH) 
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where 

(F, K{r) F) = \j ^ ^ ^) [F(^) - F(l» - r^)] \ (B12) 

Rotational invariance implies that the extremum must have the rotationally covariant form 

F s (h;p', P ,k) = hF s (h; P ',p,k). (B13) 

The number of integrations necessary to evaluate (F, K(r) F) can then be reduced by 
expanding the square in (B12), switching integration variables to the dimensionless variables 
u = h/(rmi)) and w = (h — rq L ) / '(rmc), and then performing the angular integrations 
using the explicit form (B5) for A(q ± ): 

(rm D ) 4 f°° j u 2 +w 2 (u 2 + w 2 + 1) 



(hF, Kir) hF) = T f duudww 

(4tt) 2 Jo 



l2 ~ w2 \ [((u+w) 2 + l)((u-w) 2 + l)] 1/2 
x [F(rm D u) - F(rm D w)] 2 . (B14) 

Making a finite basis expansion of the form, 

F(h) = K (m) ^ m) (h) + i F i M ^ m \h), (B15) 

m=l m=l 

an efficient set of basis elements is [30] 

(U2 /A\m-1 

$i m \h) = 1 1 ' N+2 , m=l,...,N z , (B16) 



^i (m) (fe) = ,o /Mjyi . m=l,...,M, (B17) 



where the scale A is chosen to optimize convergence as the basis set is increased. 16 Since the 
basis elements $ are only functions of h 2 /A, one can see from (B14) that matrix elements 
such as (34 m) , K(r) $[ n) ) will equal T( m-£,r) times a dimensionless function depending only 
on m, n, and z = (m^r) 2 /A. That means that one can numerically generate these functions 
just once for a given basis size (by evaluating on a fine mesh of points in z and spline 
interpolating), and then repeatedly re- use their tabulated values for calculations in different 
theories with different values of mo and A. 

For further details relating to setting up the framework for numerical evaluation, see the 
treatment of photo-emission in Ref. [30]. The solution of the integral equation is somewhat 
involved but can be performed with very good numerical accuracy (< 10~ 4 relative errors). 
Once the values of the K(r) matrix elements have been tabulated, the solution of the splitting 
functions 7^. is quite fast. Consequently, even though these splitting functions appear inside 
the integrals of the Boltzmann collision term (2.23), the numerical cost of nested integral 
equations is ultimately not very large. 



An efficient choice of A can be found by numerical search as described in Ref. [30]. However, in the present 
context, we have also found that the simple choice A — pp'rril^ + k(p' —p)ml s s works reasonably well. 
This is the value of h 2 for which the last term in Eq. (B4) equals the preceding terms. 
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APPENDIX C: ASYMPTOTIC BEHAVIOR OF INVERSE LOG EXPANSION 



Here we prove that the expansion (4.11) in inverse powers of ln(///mD) is an asymptotic 
series with vanishing radius of convergence. We also give a (non-rigorous) argument that 
the expansion coefficients have non-alternating factorial growth, implying the presence of 
singularities in the Borel transform on the positive real axis. 



1. What we need to show 



Within the i = 2 or C-odd £ = 1 symmetry channels of interest, both the linearized 
collision operator C, and its leading-log piece A defined in Eq. (4.4), are real, symmetric, 
positive definite linear operators. 17 This can be seen by looking at the explicit forms (2.21) 
and (4.5). In other words, within the symmetry channels of interest, the operators A and C 
are both invertible. 

However, we are not assured that the difference operator 5C(fi/T) which generates the 
inverse log expansion is positive definite, because its definition (4.6) involves a limiting 
procedure in which we subtract larger and larger multiples of A as we take the va^jT — > 
limit. There is no guarantee that adding a finite multiple of A to 5C(/i/T) will yield a 
positive operator, and in fact it does not. In other words, even though C is a positive 
definite operator, removing all subleading (in g) contributions by defining 

C Le leading order = C ^ + *M = ^ M<V™d) A + 5C(»/T)} , (CI) 

yields an operator which is not positive definite. To show this, we will demonstrate that 
there is a family 0( m ) of functions such that 

lim v 7 = -oo , (C2) 



o;::;i.^o;:;) 



which is enough to ensure that [ln(fi / itly)) A+5C(fi/T)] cannot be free of negative eigenvalues 
for any finite ///mo- Since A has positive spectrum, this implies that the expansion (4.10) 
in inverse powers of \n.(n/mn) must have a vanishing radius of convergence. 18 



More generally, C is positive semi-definite with zero modes associated with conserved charges: one C- 
even I = zero mode associated with energy conservation, one C-even I = 1 zero mode associated with 
momentum conservation, and various C-odd I = zero modes associated with conservation of quark 
flavors. The leading-log operator A has the same structure except for one additional I = zero mode 
associated with total particle number (gluon plus quark plus anti-quark) conservation. 
Here's a general proof. Let C(t) = A + tB, with A an Hermitian positive definite and hence invertible 
operator, B Hermitian, and t real. A necessary (but not sufficient) condition for the convergence of the 
Taylor expansion of C(i) -1 in powers of t at t = i, is the existence of C(t) _1 for all t between and t. As 
t increases from zero, C(t) first fails to be invertible when some eigenvalue crosses zero. If (y, C(t) v) < 
for some vector v, then one or more eigenvalues must have crossed zero for t < i , implying that i is outside 
the radius of convergence. 
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2. The test function sequence {0( m )} 

In the limit ttld/T — > 0, which is used to define A and SC(n/T), the leading-log contri- 
bution dominates transport coefficients. The leading-log result is associated with t-channel 
(or w-channel) 2^2 processes with momentum transfers q in the range m D C g < T. For 
the hard particle momenta p ~ k ~ T that dominate transport, one may then approximate 

Xi-j (p') = Xi-j (P+9) - Xi-j (p) + 9 • VpXi -j (p) (C3a) 

and 

Xi-j(k') = Xi-j(k-q) ~ Xi-j(k) - q ■ V k Xi-j( k ) ( C3b ) 

in the collision term (2.7). This approximation, plus the neglect of screening, reduces 
C 2 ^ 2 /(g T) to A times a logarithmically divergent integral J dq/q. One or the other of 
these simplifications break down outside the region tod C g < T, and so this integral more 
properly gives a factor of ln(T/m D ), sometimes called a Coulomb logarithm. For details, see 
Ref. [11] (or the somewhat analogous discussion of the non-relativistic case in Ref. [48]). 

Now suppose that we artificially consider functions Xi-j(p) that are similar to those one 
finds in a leading-log calculation of transport coefficients but which abruptly cut off above 
some momentum scale with m D C p, < T. Specifically, consider replacing xip) by the 
test function 

<P ( *\p) = y e- p/p * • (C4) 

The small q approximations (C3) now break down for p > rather than p > T '. There- 
fore, after making the same approximations as before, one might expect to find the same 
leading-log result for the collision matrix element except for a replacement of ln(T/mD) by 
ln(p*/wiD)- In other words, 

, ?<l$.j) = (Vft , [ln(p*/m D ) + 0(1)] . (C5) 

We will sketch a more detailed argument momentarily. Subtracting the same matrix element 
of A \n.(n/ms) from both sides, and then taking the limit m D /T — > (as dictated in the 
definition (4.6) of 5C) with p*/T held fixed, yields 

^ = ln(p,/ Ai ) + 0(l). (C6) 



Consequently, to generate a sequence of functions demonstrating (C2) one may simply take, 
for example, p* = e~ m T with m — 1, 2, 3, • • • . 



3. Details 

To justify the estimate (C5) adequately, we must show that there are no other contribu- 
tions to the matrix element which are as large as the small exchange momentum contribution 
already considered. For simplicity of presentation, we will choose to discuss explicitly the 
case where all departure from equilibrium is carried by fermions. 
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As a benchmark for further discussion, it is convenient first to note the parametric size 
of matrix elements of A. Direct evaluation of Eq. (4.5) using x(p) = 4^*\p) shows that the 
expectation value scales like p^ +1 , so that 

(*w,^w)~^ (C7) 

with comparable contributions coming from gauge boson exchange and fermion exchange 
diagrams. Now consider the contribution of 2 <-> 2 processes to the same expectation value 
of C. For definiteness, start by considering t-channel processes and work in the t-channel 
representation of (A18). The integrand of the 2^2 contribution (2.22) will be exponentially 

suppressed by the test function factor [(f>i*}.j(p) + • • • — (f)\*).j(k')] 2 unless at least one of the 
momenta p, p', k, or k' are 0(p*). This is only possible for q — \u\ < p*, and either p or 
k within of its lower integration limit. This gives an 0{p1) phase space suppression. 

For q ~ T, we will have t ~ p*T and dominantly s ~ u ~ p*T as well [because each 
Mandelstam variable can be written in the form — (P* ±P) 2 = ±p*p (1— cos#) for some one 
of the external momenta P with p ~ T\. Therefore, 19 the contribution of t-channel processes 
from the region q ~ T to the matrix element (0W,C0OO) is q^u+i j T ±y Thig ig smaller; 

by one power of p*/T, than the small-g contribution which leads to the matrix element (C7) 
for A (which also multiplies a logarithm). Contributions from smaller exchange momenta, 
g C T, will contain greater phase space suppression which can only be compensated in those 
terms whose matrix elements contain small q divergences. 

For s-channel contributions, and interference terms, one finds the same result. To see 
this explicitly, take the matrix elements of Table III with the substitution (A13) applied 
to the s 2 /tu and u 2 /st terms, so that every denominator is either a power of just s, t, or 
u (or simply constant). Terms with just s or s 2 in the denominator do not have small s 
divergences because t and u become small as well. For terms with small t or u divergences, 
we have seen that only the small q region can be non-negligible. Consider matrix elements 
that diverge as t~ n or u~ n with n = 1 or 2, and focus on the case of t. For q ~ p*, we have 
t ~ pi, which is suppressed by a single power of p* from the size t ~ p*T relevant for q ~ T. 
The contribution of these terms from q ~ p* is then enhanced by (T/p*) n ~ 1 compared to 
the q ~ T region discussed previously, with one power of p*/T reflecting the greater phase 
space suppression on q. Additionally, for the soft fermion exchange contributions to the 
qg <-> gq and qg <-> gq processes, there will be an additional Bose enhancement factor of 
order /|(p*) ~ T jp* because, for such small momentum transfers, one of the final or initial 
gluons must have 0(p*) momentum if one of the initial or final quarks does. 

For terms which do not involve soft fermion exchange, the resulting contribution to the 
matrix element of C is therefore 0(pl e ~ n+3 /T 5 ~ n ), which is (p*/T) 2 ~ n times the size (C7) 
of the matrix element of A. For qg <-> gq and q <-> gq soft fermion exchange contributions, 
there is an additional enhancement of T/p*, giving a result which is (p*/T) 1_n times (C7). 
Since the degree of divergence n is at most 2 for gluon exchange, and 1 for fermion exchange, 
this shows that even for test functions which are peaked at T, the dominant 2^2 



We can ignore the possibility of a Bose enhancement fo{p*) ~ ^Yp* from the external particle with 0(p*) 
momentum that is associated with x, because we are only considering the case where the fluctuation \ 1& 
in the fcrmionic sector. 
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matrix elements are the underlined ones of Table III — precisely those matrix elements that 
contribute to the leading-log result. The resulting contributions to {6^*\C(f)^\ are indeed 
dominated by q < p*, as claimed earlier, and the contribution from q ~ is of the same order 
as the result (C7) for A. For q <C p*, the small q divergences are softened by the cancellation 
of x factors in (2.7), which was the basis for the leading-log approximation, reviewed earlier 
in discussing (C3). As discussed previously, the contribution from mo C g C p t then 
exceeds, by a logarithm, the contribution from q ~ p*. 

Analyzing 1^2 processes is substantially more complicated, but eventually leads to the 
result that the expectation value {4>^*\C is 0{p 2i+1 /T 3 ), without any logarithmic 

dependence on p*/T or m^/T. Hence, 1 <-> 2 contributions remain sub-leading by a loga- 
rithm even when <C T. We will only outline the analysis. The factors of the test function 
(f)^ in Eq. (2.23) cause the integrand to be exponentially small unless either p or k are 
< p*. Symmetry of the integrand allows us to focus on p < p*. The dominant contribution 
comes from k ~ T (and so p' ~ T), as can be separately verified. Estimating the collision 
operator (Bll) as K s F{h) ~ g 2 m\ ) Tk 2 V 2 h F ', and noting that 5E ~ m^ ) /p+ h 2 /(pk 2 ), we 
find that the SE term dominates, so the integral equation can be solved by expanding in 
5E > IC S ~ g 2 mlk 2 T/h 2 with the result that ReF ~ K s h/(5Ef. Taking k ~ T and 
substituting into all terms in Eq. (2.23), we find that T q ~ p 2 g 4 T 5 , 7| ~ 7^ ~ (g 2 T) 2 , and 

finally obtain C 1 ^ 2 ^*)) ~ p 2e+1 /T 3 . The contribution from /c ~ is suppressed by a 

further power of p*/T. 



4. Lack of Borel summability 

Having seen that 5C is not relatively bounded by A [c.f. Eq. (C2)] due to a relative 
logarithmic enhancement at small momentum, we can obtain a (non-rigorous) estimate of 
the large order behavior of the coefficients Q n of the inverse log expansion by restricting 
the domain of 5C(fi/T) to functions whose support lies in the momentum region p < p max 
for some p max T. In other words, multiply 5C on both sides by a projection onto this 
subspace. Within this subspace, 20 

6C(fi/T) » A 1/2 [Hp/fi) + K] A 1/2 , (C8) 



In more detail, let (p^ 1 ^ represent any function roughly similar to (C4) with ~ pi — that is, a function 
which is cut off for momenta large compared to p± <C T. Now consider matrix elements of SC of the 
form (<j)( pi \ SC </>( p2 )). One can repeat the analysis of Appendix C2 with the only new ingredient being 
that the matrix element need not be diagonal. However SC is diagonally-dominant: if pi <p 2 , then the 
normalized off-diagonal matrix element ((t> [pi \ (p2) )/[(0 (pi) , </> (pi) ) (4> {p2 \ (j) (p2} )] 1/2 is small compared 
to the normalized diagonal elements {4>^ Pi \ SC (j)^^) / ((j)( pi \ 0^ Pi ') for i = 1,2 [as can be seen from the 
explicit form (4.5)]. When p\ ~ P2, the momenta p that dominate the integrals representing the matrix 
element will be p ~ p\ ~ p2- The logarithm factor in the discussion of Appendix C2 will become 
ln(pi///) ~ ln(p2/ 'a 4 ) and w e have simply replaced this by ln(p/(i), given which momenta p dominate. The 
particular ordering of factors in the expression (C8) is not important at this level of precision; the given 
form is Hermitian (as it must be) and convenient for what follows. 
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for some constant K which depends on just how the cutoff is implemented. Inserting this 
approximation into the result (4.12) for the expansion coefficients Q n , one has 

Q n+ M « Q^r(») = \ (-i) n (Si-j, A- 1 ' 2 V(p) [Hp/») + K] n A^ 2 S t ..,) , (C9) 

where V(p) = 6(p m ax — p) denotes the projection operator restricting the domain of 5C. 
Assuming that p max < e~ K /x, the approximation (C8) to 5C is negative-definite, implying 
that the expression (C9) for Q n+ i(n) is strictly positive. Given this positivity, we can apply 
the same variational approach discussed in section II D. Specifically, using the test functions 
4>(*\ one has the lower bound 

USi-j, *&) (o^- 5,.,) 

QnZriti > -T-, K — A - - \ M , ■ (CIO) 

(0« , aw [- in( P //i) - k] n A^ 

Consider now, for simplicity, the case of I = 1 with all departure from equilibrium carried 
by fermions, as is relevant for flavor diffusion. An easy estimate shows that 

(S,^) ~pj, (Cll) 

while from Eq. (C7) one has 

^t 1 / 2 [- ln(p//i) - K] - n A 1 ' 2 ~ p2 [- ln(^//i) - #i ] ~ n , (C12) 
for some new constant K\. Hence 

QTr^) > K2 pi [- ln(p,//i) - K,] n , (C13) 

for some positive constant i^ 2 - (The 0(1) constants K x and K 2 do not depend on p* or ra.) 
Since this is a lower bound, we may maximize over all values of p* (subject to p* < e~ Kl //), 
which gives p*| op t = e _ri//5 e~ Kl /i and results in the estimate 

Q— (/x) > K 2 (e"%) 5 e"" . (C14) 

The case of shear viscosity, where £ = 2 and bosons are also out of equilibrium, is similar. 
The non-alternating factorial growth of (C14) with n immediately implies that the inverse 
log expansion (4.11) has zero radius of convergence, and moreover that its Borel transform 
has a singularity on the positive real axis (at 5). Although the estimate (C14) is based on 
the approximation (C8), we do not see how inclusion of the non-logarithmically enhanced 
parts of 5C can affect the n\ growth of the above estimate. Assuming that this estimate does 
give the correct large order behavior of the inverse log expansion, it is worth noting that 
the resulting renormalon ambiguity associated with the singularity of the Borel transform 
would only be of order (mv/T) 5 ~ g 5 . 
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